A large-A/^ analysis of the local quantum critical point and the spin-liquid phase 
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We study analytically the Kondo lattice model with an additional nearest-neighbor antiferromag- 
netic interaction in the framework of large-A/" theory. We find that there is a local quantum critical 
point between two phases, a normal Fermi-liquid and a spin-liquid in which the spins are decou- 
pled from the conduction electrons. The local spin susceptibility displays a power-law divergence 
throughout the spin liquid phase. We check the reliability of the large-A/" results by solving by 
quantum Monte Carlo simulation the A/" = 2 spin-liquid problem with no conduction electrons and 
find qualitative agreement. We show that the spin-liquid phase is unstable at low temperatures, 
suggestive of a first-order transition to an ordered phase. 
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The effects of the vicinity to a quantum critical point 
(QCP) on the low-temperature physical properties of var- 
ious systems are presently attracting great attention [1]. 
There is in particular growing consensus that the anoma- 
lous metallic properties observed in some heavy-fermion 
materials as well as the presence of superconductivity 
near magnetic phases could be interpreted in terms of 
the proximity to a QCP [2]. The issue of QCPs in super- 
conducting cuprates is also a hotly debated subject [3]. 

In the traditional picture of the QCP in metallic sys- 
tems [4,5] collective spin-fluctuations of specific wavevec- 
tors become critical at the QCP. A substantially differ- 
ent scenario has recently been suggested [6] in which the 
quantum phase transition involves all wavevectors simul- 
taneously. At this local QCP the local degrees of freedom 
develop long-range correlations in time [8] leading to a di- 
vergent local susceptibility. At the onset of magnetism 
the energy scale that characterizes the coherent metallic 
phase vanishes. 

It was proposed [6] that this peculiar type of behavior 
arises in a particular Kondo-lattice model with nearest- 
neighbor magnetic RKKY interactions. The model was 
studied within Extended Dynamical Mean-Field Theory 
(EDMFT) which reduces the problem to an effective im- 
purity model in which a localized spin interacts with 
fermionic and bosonic baths [10,11]. However, some im- 
portant aspects of the self-consistent procedure used in 
Ref. [6] to close the set of EDMFT equations require fur- 
ther scrutiny [8]. 

In this paper we use a large-A/" approach to solve 
analytically the Kondo-lattice-RKKY model of Ref. [6] 
fully implementing the self-consistent EDMFT scheme. 
We supplement this large-A/" analysis with a Quan- 
tum Monte Carlo (QMC) investigation of the Kondo- 
decoupled phase of interacting spins represented by the 
pure RKKY model. 

We start from the following Hamiltonian, 
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The operators c\ a and Si represent respectively the 
creation of a conduction electron of spin a =|,J, and 
a localized spin at the zth site of a regular lattice of 
size N. We denote by Sj the local spin density of 
the conducting electrons. The hopping terms Uj corre- 
spond to a non interacting electronic density of states 
p°(e) = 1/N^ fc <5(e — e*) and the antiferromagnetic 
RKKY couplings are described by a spectral density 
pi (e) = 1/N J2 q 3{ e ~Jq) where eu and J q are the Fourier 
transforms of the nearest neighbors couplings t\j and Jij , 
respectively. We assume for simplicity that p° c and pi are 
even functions that vanish outside the intervals [-D, D] 
and [—J, J], respectively. 

In the EDMFT approach the above model is imple- 
mented on a lattice with a large coordination number z, 
so that both the hopping and the magnetic RKKY cou- 
pling are scaled as Uj — > Uj / s/z, Jij — » Jij / \fz. 

To leading order in a large- z expansion we obtain the 
following effective action for the local degrees of freedom 
on the 0-th site (say) [13] : 
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where the correlators Qq 1 (r) and x^ 1 (t) are self- 
consistently determined cavity fields [9-11]. 

There are two steps in the implementation of the 
EDMFT procedure [9-11]. The first (and usually most 
difficult) one consists in solving the local impurity prob- 
lem for fixed Qq and x^" 1 in order to find the spin- 
symmetric local electronic Green function G c (t) = 



1 



-(Tc ff (r)4 (0)) and the local magnetic susceptibility 
XIoc(t) = (TS(r) • S(0)). Next, the condition that the 
impurity site is equivalent to any other lattice site is im- 
posed through a set of self-consistency relations that al- 
low us to express the cavity fields Q^ 1 and Xq 1 m terms 
of G c and xioc- 

Here we implement this procedure in the framework 
of a large- TV expansion [12]. The SU(2) spins arc re- 
placed by SU(7V) operators in the fermionic representa- 
tion S a ,a' = flU> - <W/2 with a, a' = 1,...,7V. The 
/-fcrmions are subject to the constraint X^Li fafa = 
Af/2, enforced through the introduction of a Lagrange- 
multiplier i\(r) [12]. The coupling constants must now 
be rescaled so that physical quantities remain finite in the 
TV — > oo limit. A consistent rescaling of all the couplings 
in the Hamiltonian (1) is not possible without loosing the 
dynamical character of the spin fluctuations. We adopt 
instead a different approach in which we rescale the cou- 
plings in the action (2) as Jk — > Jk/N and Xq" 1 — > 
Xo 1 /Af, and define the local susceptibility per spin com- 

poncnt Xloc (r) = Af~ 2 E„*> (fa(r)fa> (r)/J, (0)^(0)). 
With this rescaling all the terms in Eq. (2) are 0{Af). 

In the J\f — > oo limit, the hybridization parameter 
r = —Jk{cI f a ), the Lagrange multiplier A and the chem- 
ical potential p are determined from the saddle point 
conditions that may be written in the compact form 

~,^y} = -{G/c,G/,G c }(t = 0-). (3) 

G f (r) = -(T/ CT (r)/t(0)), G c (r) = -(Tc CT (r) 4(0)) and 
G/ C (r) = — (T/ CT (r)cJ.(0)) are the dressed /-fermion, con- 
duction electron and mixed Green functions. The ef- 
fect of spin fluctuations is embodied in a self-energy cor- 
rection £/(t) = — 2xo 1 (t)G/(— t) to the bare / prop- 
agator, t/j 1 (o>„) = iuj n + X — £/(u;„). The dressed 
Green functions are G c (u> n ) = G a c (iu) n + p, — r 2 Qf(co n )), 
Gf(u) n ) = Qf{oo n ) + r 2 Gj(u n )G c (u) n ) and G fc (u> n ) = 
rQf{<jj n )G c {w n ), respectively. The usually difficult step 
of obtaining xioc from the bath Xo 1 can easily be per- 
formed in the large- Af limit in which 



Xioc(r) = -G f (r)G f (-r). 
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The parameter r describes the binding of the c elec- 
trons to the / ("spin") degrees of freedom. If r is fi- 
nite there is Kondo compensation of the localized spin. 
For sufficiently small values of J we find a Fermi-liquid 
(FL) phase with r ^ for temperatures T < Tk{J), 
the Kondo temperature of the system. In this low-T and 
low- J region, the FL is characterized by a finite coher- 
ence energy scale cfl cx Tk and a large Fermi surface 
containing 1 + n c states (Luttinger theorem) . The phys- 
ical properties of this heavy Fermi liquid are similar to 
those discussed in Ref. [13]. In particular, the local mag- 
netic susceptibility xioc(T = 0,u = 0) remains finite. 



Upon increasing J, the physical quantities characteriz- 
ing the Kondo-screened FL (r, Tk, and epi) gradually 
decrease and vanish at a critical value of the magnetic 
coupling J c <~ T K , the Kondo temperature for J = 0. At 
the same time xioc{T = 0, u> = 0) increases continuously 
and diverges at J = J c . The paramagnetic solution with 
r = can be described in terms of a gas of free c electrons 
with a small Fermi surface enclosing n c states decoupled 
from a spin-liquid (SL) of strongly correlated / fermions. 
Notice that, since our large- Af approach neglects the fluc- 
tuations of r, we cannot distinguish between the critical 
behavior at the QCP and that inside the SL phase. 

Of course, a global magnetic instability may occur at 
a finite temperature T — T c before this local instability 
of the Fermi surface has a chance to develop. This will 
happen if x(<7 = Q,w = 0, T) diverges at T = T c , where 
Jq = —J corresponds to the lower edge of the RKKY 
spectral density pj. In this case, the local QCP behavior 
will be hidden by the magnetic instability but its effects 
may still be observable provided that spin ordering is 
established at sufficiently low temperatures. 

In order to proceed we must fix the actual prescription 
used to close the self-consistency loop relating Xo 1 to 
Xioc- Following the EDMFT scheme of Refs. [6,10,11], 
we introduce the local conduction electron and mag- 
netic self-energies, S c (w„) = tJ _1 (w n ) — G ( 7 1 (w„), and 
M(uj n ) = Xo : K) + Xi c{u n ). These quantities may be 
expressed in terms of O(J\f ) physical quantities as 
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with x 1 (q,u„) = M(uj n ) + J q and G c 1 (k,uj n ) = iv n - 

S c (w„) - Cfc. 

The results depend crucially on the behavior of pi{e) 
near the band edge. We discuss first the case in which the 
RKKY spectral density pj (e) is finite at the lower edge of 
the band, p(—J) = pi, a situation that arises when the 
spin correlations have a two-dimensional character [6,10]. 
Then we may explicitly integrate Eq. (6) to find 



Xloc(u n ) = P/ln 



M(to r , 



J 



M{uj n ) - J 
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where lZ(oj n ) is a regular function of frequency. Notice 
that in this case a second-order magnetic transition [sig- 
naled by M(u) n = 0) = J] is necessarily accompanied by 
a divergence of xioc- Such a divergence at T = can 
only take place when r vanishes, that is at the quantum 
critical point, J = J c , and inside the SL phase J > J c . 

A complete solution of the model at the QCP and in 
the SL phase can be obtained in the large- Af approach 
by solving Eqs. (3) and (4) for r — supplemented with 
the relation 
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which follows easily from Eq. (7). 

Since the solution obtained from the RG treatment 
for the local impurity problem [6] yields a power-law 
behavior for M and a logarithmic divergence for \ioc 
at low u> and T, we first check the validity of these 
results within our large- Af scheme. We thus assume 
that, at low frequency, M(w n ) — J ~ (|o;„|/J) 7 and 
Xioci^n) = 7 At hi ( J/\u> n \). Analytic continuation to real 
frequencies leads to: 



Im[ Xo -V)] = 



7rsgn(cj) 
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The imaginary part of the self-energy for real fre- 

quencies is 



/-t-oc 
delm [xo 1 ^)} p f (u-e) x 
-OO 

[n B {-e)+n F {oj-e)} . (10) 

In the low-frequency regime, with the (self-consistently 
verified) assumption that dominates lo, the 

corresponding fermionic Green function is given by 
Gf(u) w -EJ 1 ^). Taking then p f (uj) = 
^/"fpi/ (2tt) In I J/u)\/y/\co\ [corresponding to £y(u>) = 
— -y/jcjf/ (y/jpi hi\J/u>\)], and with the correlator of 
Eq. (9), one can check that Eq. (10) is verified. 
Using these results in Eq. (4) we find xioci^n) — 
(7P//67r) [In | J/uj 71 \] 3 . We now notice that this finding 
is inconsistent with the previously assumed simply log- 
arithmic behavior of xioc- On the other hand, one can 
further check that an Ansatz pf(u>) ~ 1/ y/\Ju)\ (neces- 
sary to get a purely logarithmic xioc) does not verify 
Eq. (10). All the above shows that, within our large- Af 
scheme, a correlator of the form x " 1 (o;, i ) oc LrG (J/w n ) 
produces additional logarithmic corrections to the local 
susceptibility, which prevent self-consistency. 

On the other hand, one can still find self-consistent 
solutions of the problem. Starting from the Ansatz 
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XtocK) = {A/ J) [\w n \/J] 
G f (u n ) = -i(B/J) [K|/J] 



(11) 
(12) 



diverging at low frequencies, one finds ImXo [^) ~ 
( J/A)(\u)\/ J) s sgn(uj) sin(<$7r/2). Inserting these expres- 
sions in Eq. (10) one can determine (T — 0,u > 0) 



£';(w) = -(2JB/nA) (uj/J) 5 ' a+1 

X sin (tt6/2) cos (ira/2) f dxx s /(l - x) a . (13) 
Jo 

Assuming that 5 < a (that will be self-consistently veri- 
fied), one has Gf(u>) w — S7 (a;). Compatibility between 
Eqs. (12) and (13) requires 
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FIG. 1. Quantum Monte Carlo results for the decoupled 
model. We used up to 512 Trotter time-slices and 1.5 x 10 5 MC 
steps/time-slice, (a): Jxiocij) for J/T = 50. Ao ~ 0.02059 is 
a short-time cutoff obtained fitting the data to the functional 
form shown in the legend, (b): T-dependence of the static 
local susceptibility 
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On the other hand, self-consistency requires that xioc 
calculated from Eq. (4) for the specific G/ considered 
here coincides with Eq. (11). This leads to the additional 
condition 
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where the last equality follows from Eq. (15). Surpris- 
ingly, this equation (in conjunction with Eq. (14)) can be 
solved analytically. We find 5 = 2a — 1 = 1/3. We have 
tested these large- Af findings by a numerical approach 
based on the QMC algorithm of Ref. [14]. Since no algo- 
rithm exists capable of treating the Kondo and RKKY 
interactions on the same footing, we confine ourselves to 
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the investigation of the paramagnetic Kondo-decouplcd 
SL phase. Therefore we solve the self-consistency equa- 
tions for the pure RKKY model represented by the last 
term of the action Eq. (2) for standard SU(2) spins with 
a constant pi(e) = 1/(2 J). The local susceptibility in 
imaginary time thus obtained is shown in Fig. 1(a). The 
low-T, long-time behavior deduced from the numerical 
results is XIoc(t) ~ [sin(7rr//3)] _1 / 3 . This corresponds to 
the following scaling form 

X'LM ~ J _1 (w/J)" 5 ^(w/T) (18) 

with T s {x) = x s \T(^ +^)| 2 sinh(f) and S = 2/3. 
Eq. (18) gives a power-law divergence at low energy, 
X'UT = 0,w) ~ J-iuj-i and x[ oc (T,u; = 0) ~ 
J~zT~3 . The T-dcpcndcnce of x/ oc (w = 0), displayed 
in Fig. 1(b), shows indeed this behavior. Quite interest- 
ingly, apart from the difference in the value of the expo- 
nent S (that was expected as critical exponents do depend 
on AT) the "exact" QMC results for SU(2) spins are quali- 
tatively similar to those of the simpler large- AT approach. 
We expect on general grounds that the inclusions of 1/Af 
corrections should improve the agreement between the 
1/Af expansion and QMC results. It is important to no- 
tice that if the solution (18) for \ioc is used to determine 
M(uj) from Eq. (8), the resulting M"(u) exhibits oscil- 
lations below a temperature T*, a behavior that violates 
the condition of thermodynamic stability uM"(ui) < 0. 
We have determined T* in two different ways, from our 
analytic expressions in Eqs. (8) and (18), and by analyt- 
ically continuing our numerical results on the imaginary 
frequency axis using the method of Pade approximants. 
The two methods yield T* ss 0.03J. This is proportional 
to, but substantially smaller than J, so that a wide range 
T* < T « J of temperatures still exists within which 
the power-law given by Eq.(18) holds. The instability 
found at T* could signal that the magnetic transition is 
first-order as suggested by other authors [15] (see also 
Rcf. [16]); we will address this issue elsewhere [17]. 

We discuss next the case in which pi(e) has a square 
root singularity at the band edge which corresponds to a 
three dimensional spin-wave spectral density. A partic- 
ularly simple case is that of the semi-circular density of 
states, pi(e) = (Jir/2)~ 1 sjl — (e/J) 2 . Then it is easy to 
show that Xq 1 (w„) = (J 2 /4)x; c(w„) [9]. This particular 
model was studied in Ref. [13] in the context of the SL 
state first discussed in Ref. [18]. At the local QCP, the 
fcrmionic Green function diverges like Gf(u>) ~ 1/ \J J|w| 
which yields x*oc(w) ~ J -1 ln( J/\uj\) allowing to close 
the self-consistency procedure without the appearance of 
additional logarithmic corrections. These results hold 
as long as the magnetic self-energy M(iv n ) remains out- 
side the interval [—J, J], so that no magnetic instabil- 
ity occurs. However, in this case M{ui n = 0) = J for 
Jxiod^n = 0) = 2, a condition satisfied at a finite tem- 
perature T c <~ J. Since at the transition J < T^-, it oc- 



curs while Tk (and r) are still finite. Then, as pointed out 
in Ref. [6], in the three dimensional case the local QCP 
is masked by the finite temperature magnetic transition. 
The analysis of this QCP, expected to be of the Millis- 
Hertz type [4,5], lies beyond the power of our lowest-order 
large- N theory 

In summary, we carried out a largc-A/" implementa- 
tion of the EDMFT equations for the Kondo-lattice- 
RKKY model. This simple scheme allows us to find self- 
consistent solutions of the EDMFT equations without 
any ad hoc input. We find in agreement with Ref. [6] 
that there exists a local QCP separating a FL and a SL 
phase. However, this approximation, which neglects the 
fluctuations of r, does not allow us to distinguish be- 
tween the properties of the QCP and those of the SL 
phase which correspond respectively to the unstable and 
stable fixed points of the RG equations of Ref. [6] . We 
showed that the large- Af results for the SL phase are reli- 
able by comparing them with those of an exact numerical 
solution of the SL problem for Af = 2. Both our meth- 
ods (in agreement with the RG results of Ref. [7]) show 
that a correlator Xq" 1 <~ ui s corresponds to a local sus- 
ceptibility xioc ~ uo~ s . On the other hand we find that a 
correlator x^ 1 ~ m_1 w does not correspond to a simply 
logarithmic susceptibility, showing that the limit 6 — > 
can not be taken trivially. Although this conclusion can 
only be safely drawn in connection with the SL problem, 
it suggests that caution must be exerted when taking the 
same limit at the QCP [19]. 

Finally, we showed that the power-law solution that 
we found in the SL phase is thermodynamically unsta- 
ble below a temperature T* <C J- This suggests that a 
first-order phase transition could take place at low tem- 
perature as suggested by others [15]. 
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